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ABSTRACT 

Anomalous X-ray pulsars, compact non-pulsing X-ray sources in supernova remnants, and X-ray 
bursters are three distinct types of sources for which there are viable models that attribute their X- 
ray emission to thermal emission from the surface of a neutron star. Inferring the surface area of the 
emitting regions in such systems is crucial in assessing the viability of different models and in providing 
bounds on the radii of neutron stars. We show that the spectroscopically inferred areas of the emitting 
regions may be over- or under- estimated by a factor of <^ 2, because of the three-dimensional geometry 
of the system and general relativistic light deflection, combined with the effects of phase averaging. Such 
effects make the determination of neutron-star radii uncertain, especially when compared to the ~ 5% 
level required for constraining the equation of state of neutron-star matter. We also note that, for a 
given spectral shape, the inferred source luminosities and pulse fractions are anticorrelated because they 
depend on the same properties of the emitting regions, namely their sizes and orientations. As a result, 
brighter sources have on average weaker pulsation amplitudes than fainter sources. We argue that this 
property can be used as a diagnostic tool in distinguishing between different spectral models. As an 
example, we show that the high inferred pulse fraction and brightness of the pulsar RXS J1708— 40 are 
inconsistent with isotropic thermal emission from a neutron-star surface. Finally, we discuss the implica- 
tion of our results for surveys in the soft X-rays for young, cooling neutron stars in supernova remnants 
and show that the absence of detectable pulsations from the compact source at the center of Cas A (at 
a level of ^ 30%) is not a strong argument against its identification with a spinning neutron star. 

Subject headings: relativity — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Neutron stars in nature appear in various different fla- 
vors. They were first unambiguously discovered in radio 
wavelengths as rotation-powered pulsars and later on, in 
X-rays, as accretion-powered pulsars and bursters. Such 
systems, however, constitute only a small fraction of the 
expected number of neutron stars in the galaxy, as in- 
ferred from estimates of supernova rates (e.g., Kaspi 2000). 
As a result, ongoing searches exist for detecting neutron 
stars in different manifestations, e.g., as isolated stars ac- 
creting from the interstellar medium (e.g., Belloni et al. 
1997) or simply cooling by thermal emission (e.g., Pavlov 
et al. 1996). Different emission mechanisms are thought 
to operate in different types of neutron stars and possibly 
even between different wavelength bands in a given sys- 
tem. For at least three distinct classes of neutron stars, 
which we discuss below, viable models exist in which radi- 
ation emerging directly from their surfaces is responsible 
for their high-energy spectra. 

A number of slow X-ray pulsars, often called the anoma- 
lous X-ray pulsars (AXPs; Mereghetti & Stella 1995) form 
a distinct class of neutron stars with soft X-ray spectra 
but no radio or optical counterparts. Because of their high 
spin-down rates and association with supernova remnants 
(SNRs) they are thought to be young and related to the 
soft 7-ray repeaters, in their quiescent states (see Hurley 
2000). It is still an open question whether such objects are 
neutron stars accreting from a fossil disk (e.g., van Paradijs 
et al. 1996; Chaterjee, Narayan, & Hernquist 2000), are 



powered by magnetospheric emission (e.g., Thompson & 
Duncan 1995, 1996), or even emitting thermally (Heyl & 
Hernquist 1998). 

Another class of objects, potentially related to the AXPs 
are the compact, non-pulsing, soft X-ray sources that are 
being discovered within SNRs. The most recent of such 
compact objects is the central source in Cas A discovered 
with the Chandra X-ray Observatory (Tananbaum 1999; 
see also Pavlov et al. 2000; Chakrabarty et al. 2000). Their 
spectra and luminosities are typical of what is expected for 
cooling, young neutron stars (see also Pavlov et al. 1996; 
Zavlin et al. 1998, 1999). However, the absence of co- 
herent pulsations in their X-ray brightness allows for the 
possibility that they are black holes accreting from fall- 
back material (see, e.g., Pavlov et al. 2000; Chakrabarty 
et al. 2000). 

Finally, weakly-magnetic (< 10 10 G) accreting neu- 
tron stars often show thermal emission from their sur- 
faces during thermonuclear flashes, the so-called Type I 
X-ray bursts (see Lewin, van Paradijs, & Taam 1996). In 
six bursters, relatively coherent oscillations are detected 
during the bursts, at frequencies ~ 300 Hz (see, e.g., 
Strohmayer et al. 1996). The coherence and stability of 
the oscillations suggests that they occur at the neutron- 
star spin frequencies and are caused by the non-uniform 
pattern of burning in their surface layers. 

Understanding the properties of the X-ray emission from 
such thermally emitting neutron stars, especially in con- 
nection to the presence or absence of detectable coherent 
pulsations, is crucial in assessing the viability of different 
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models. Moreover, comparing model spectra to observa- 
tions offers the possibility of measuring the radii of neutron 
stars and hence constraining the properties of neutron-star 
matter (see Lattimer & Prakash 2000 for a recent discus- 
sion). Recently, there has been significant progress in cal- 
culating spectra emerging from neutron-star atmospheres 
with various compositions (see, e.g., Pavlov et al. 1996; 
Rajagopal, Romani, & Miller 1997). Constraints on the 
amplitudes of oscillations from compact stars with non- 
uniform surface emission have also been studied in con- 
nection to the observed oscillation amplitudes in cooling 
radio pulsars (Page 1995), bursters (see, e.g., Miller, & 
Lamb 1998; Weinberg, Miller, & Lamb 2000), and AXPs 
(DeDeo, Psaltis, & Narayan 2000). 

In this article, we address a number of issues related 
to the emission from a spinning neutron star with a non- 
uniform surface brightness. In particular, in §2 we study 
the effects of the three-dimensional geometry of the sys- 
tems, phase averaging, and the general relativistic deflec- 
tion of light. In §3, we demonstrate that, when the ra- 
diation emerges from a localized surface area on a rotat- 
ing star, the radiation flux reaching an observer at infin- 
ity may be substantially different compared to the case of 
isotropic, spherically symmetric emission from a compact 
star with the same surface area (see also Zavlin, Shibanov, 
& Pavlov 1995; Leahy & Li 1995). As a result, when such 
effects are not taken explicitly into account, the hot-spot 
sizes, inferred from the observed fluxes and temperatures, 
may be significantly over- or under-estimated. We show 
that general-relativistic deflection typically reduces this 
discrepancy, depending on the compactness of the neutron 
star. 

We also note that the surface area on the neutron star, 
from which the localized emission emerges, determines 
both the flux that reaches the observer at infinity and 
the amplitude of pulsations at the stellar spin frequency. 
In fact, for stars with the same local surface temperature 
and emerging spectrum, an anticorrelation is expected be- 
tween the luminosity and pulsation amplitude. Therefore, 
simultaneous consideration of these two properties for a 
given system can offer strong constraints on the operating 
emission mechanisms. In §4, we investigate this property 
and its implications for the observations of compact X-ray 
sources with and without detectable pulsations. 

2. THE BRIGHTNESS OF A SPINNING NEUTRON STAR 

In this section, we calculate the phase-averaged radi- 
ation flux that reaches an observer at infinity from the 
surface of a spinning neutron star (or other compact ob- 
ject) with a non-uniform surface brightness, following the 
procedure outlined by Pechenick, Ftaclas, & Cohen (1983; 
see DeDeo et al. 2000 for the details of our implementa- 
tion). Throughout this article, we neglect any interaction 
of photons with matter between the surface of the neutron 
star and the observer, as well as any polarization effects 
(see, e.g., Shaviv, Heyl, & Lithwick 1999). We also as- 
sume that the star is slowly rotating, so that its spacetime 
is described by the Schwarzschild metric. Finally, we set 
c = G = 1, where c is the speed of light and G is the 
gravitational constant. 

We specify, as a boundary condition, the specific inten- 
sity 1(9, cj>, 9') of radiation emerging from the neutron-star 
surface, evaluated at the local rest frame and integrated 



over all photon energies. The polar coordinates (9, 4>) de- 
termine the position on the stellar surface with respect 
to the rotation axis, while the angle 9' is measured lo- 
cally with respect to the radial direction. Hereafter, we 
assume that the emission is isotropic and hence that the 
specific intensity is independent of 9'. This choice leads to 
the strongest effects of general relativistic light deflection, 
even though it is not necessarily the most appropriate for 
thermal emission from a neutron star (see, e.g., Zavlin et 
al. 1998). Calculating in detail the beaming of radiation 
requires the knowledge of the temperature stratification of 
the neutron-star atmosphere and the solution of the result- 
ing radiative-transfer problem, which is beyond the scope 
of this paper but will reported elsewhere. 

The brightness distribution on the neutron star surface 
depends on the emission process under consideration (i.e., 
thermal cooling versus localized nuclear burning) , the sur- 
face profile of its magnetic field (see, e.g., Heyl & Hcrnquist 
1998), and the possible presence of lateral metallicity gra- 
dients (see, e.g., Pavlov et al. 2000). For the purposes of 
our analysis, considering a geometry of two antipodal hot 
spots of variable size with uniform surface brightness cap- 
tures the relevant properties of these different processes. 
Following the notation used for radio pulsars (e.g., Lyne 
& Graham-Smith 1990), we denote the half-opening angle 
of each spot by p, its angular distance from the rotation 
pole by a, and the emerging constant specific intensity by 
Ins- 

The flux measured by an observer at distance d, whose 
polar coordinates with respect to the stellar rotation axis 
are denoted by ((3, $), is given by (Pechenick et al. 1983, 
eq. [3.15]) 

x (V-Soo) / h[9(x,f3);p,9 ]xdx . (1) 
Jo 

Here, Mns and i?NS are the neutron-star mass and radius, 
<7oo = — (1 — 2Mns/-Rns) f° r a Schwarzschild spacetime, 
imax = (-Rns/-^ns)/V 500, and the function h(9;p,9 ) 
is defined in Pechenick et al. (1983, eq. [3.3]). The angle 
9q measures the distance on the stellar surface between 
the center of one of the hot spots and the direction to the 
observer and is given by 

cos 9o = sin a sin (3 cos 3> + cos a cos (3 . (2) 

Equation (|l|) is simply the integral of specific intensities, 
at the distance of the observer, over the impact parameters 
b = iMns of rays that are parallel at radial infinity. The 
angle 9{x, (3) at which each parallel ray intersects the stel- 
lar surface is found by integrating the photon trajectory 
from radial infinity to the stellar surface, i.e., 

9= / [x~ 2 -{l-2u)u 2 ]- 1,2 du (3) 

Jo 

(Pechenick et al. 1983). 

The average flux, measured by an observer at infinity, 
averaged over a time longer than the rotational period of 
the neutron star, simply is 

Foo = tt/ns (V-Soo) A(a,P,p,p) ( j , (4) 
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where p = J?ns/2Mns and we have defined 



A(a,/3,p,p) 
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ft.[(9(a;, (3); a, 9n]xdx 
(5) 



When the local radiation spectrum emerging from the stel- 
lar surface is that of a blackbody of temperature Tns , then 
t^ns = f ^NS' where CT is the Stefan-Boltzmann constant. 

If the emission were spherically symmetric, an ob- 
server at radial infinity would measure a specific inten- 
sity Ioc = (— goo) Ins and a blackbody temperature = 
(V~5oo)^NS- In the absence of any prior information 
about the opening angle of the hot spots and their ori- 
entation, an observer at infinity would therefore infer for 
the emitting region a surface area 



And' 



oTt. 



(6) 



It is then customary to assume a given neutron-star mass 
and radius and correct for the effect of gravitational red- 
shifts as (see, e.g., Lattimer & Prakash 2000) 



Sinf = (-goo)Soc = (-.9oo)47T(i 2 



(7) 



Given that the real surface area of two polar caps of half- 
opening angle p is S pc = Aw(l — cos p)R^ s , we conclude 
that the error in the estimate of the emitting area is 



inf 



1 



pc 



cos p 



-A(a,f3,p,p) 



(8) 



As a consequence, estimating the radius of each of the two 



polar caps according to 2nD^ 



Si n {, leads to an error 



with respect to the real polar cap radius D pc 
to 



Anf 



l2A(a,0,p,p) 



pi?Ns equal 



(9) 



Finally, assuming that the whole neutron star is emitting 
uniformly and using equation (j7|) to infer its radius, results 
in an error equal to 



inf 



A(a,f3,p,p) . 



(10) 



3. RESULTS 



3.1. Uncertainties in the Inferred Surface Areas of the 
Emission Regions 

We present in this section the uncertainties introduced 
to the estimates of the spectroscopically inferred polar- 
cap surface areas by the three-dimensional geometry of 
the problem, the effects of phase averaging, and the gen- 
eral relativistic light deflection. 

Figures 1 and 2 show the ratio of the inferred to the in- 
trinsic surface areas of the polar caps, for different opening 
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Fig. 1. — Ratio of spectroscopically inferred polar-cap surface 
area 5i n f to their intrinsic a function of the orientation 

of the polar caps (a), for various values of the observer inclination 
(/3) with respect to the rotational axis and for different neutron-star 
radii. The opening angle p of each polar cap is held fixed at 10 
degrees. 
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2. — Same as in Figure 1 but as a function of the opening 
p of each polar cap. 



angles, neutron-star radii, and orientations of the rotation 
axis with respect to the polar caps and the observer. For 
neutron-stars that are not very relativistic (p — 4), the 
inferred surface areas can be significantly over- or under- 
estimated, depending on the relative orientations of the 
polar caps and the observer. In order to understand this 
effect, we calculate explicitly the ratio Si n {/S pc for the 
limiting case of a Newtonian star, an infinitesimally small 
emitting area, and two specific orientations. 

When a = (3 = 0°, one polar cap always appears at the 
geometric center of the stellar disk, and the flux measured 
at infinity is time independent. In this case, cos(9o = 1 
and, by symmetry, 

M0(x,/3);p^ O,0o = O°] =2vr. (11) 

Using these values in evaluating the integral (||), we obtain 

A(a = /3 = 0°,p^0,p^oo) = 2(1- cos p) (12) 
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which gives Si n {/S pc — 2. When, on the other hand, 
a = 0° and (3 = 90°, cos6* = 0, the polar caps are being 
viewed only at grazing angles. In this case, 



fc[0(x,/?);p->O,0o = 9O°]x^O, (13) 



and 



A{a = 0°,f3 = 90°,p^0,p^ oo) (14) 

and hence SW/Spc — > 0. As a result, for a Newtonian star 
and an infinitesimally small emitting area, 



< 5^ < v . 



5, 



(15) 



pc 



For increasingly more compact neutron stars, i.e., for 
decreasing values of p, the error in the estimate of the 
emitting area decreases substantially, reaching <J 5% when 
p = 2. This is the result of the strong gravitational light 
bending near the neutron star surface, which efficiently re- 
distributes the emitting photons to almost all directions of 
propagation, mimicking the spherically symmetric case. 

Figure 3 shows the maximum and minimum of the ra- 
tio S'inf/S'pc for all possible orientations of the hot spot 
and the observer, as a function of the fractional emitting 
surface area (S , pc /47ri?NS = 2f2 pc /47r) of the neutron-star 
surface, for different neutron-star radii. As expected, in all 
cases, the maximum corresponds to a — (3 = 0° and the 
minimum to a = 0° and (3 = 90°, and they both converge 
to unity as Q pc — > 2ir. However, for small polar caps and 
typical neutron-star masses and radii, there is a factor of 
^ 2 spread in the systematic uncertainty in the estimated 
area of the emitting region. 

3.2. Trends Between the Pulse Fractions and 
Luminosities of Thermally Emitting Neutron Stars 

The X-ray brightness of a spinning neutron-star with a 
non-uniform surface emission shows pulsations at the stel- 
lar spin frequency and its harmonics. Just as in the case of 
the flux observed at infinity, the amplitude of pulsations 
as well as the harmonic structure also depend strongly on 
the brightness distribution on the stellar surface and the 
degree of gravitational light bending. Therefore, inevitable 
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Fig. 5. — Maximum allowed pulse fraction from a spin- 
ning neutron star as a function of its apparent brightness [I = 
(47rci 2 F oo /10 36 erg s _1 )/(Too/0.5 keV)] inferred by an observer at 
infinity, assuming that the emission is spherically symmetric and 
without applying any redshift corrections. The local radiation spec- 
trum is assumed to be that of a blackbody and the radius of the 
neutron star is fixed at 10 km. Different relativity parameters corre- 
spond to different neutron-star masses. The data points correspond 
to the X-ray pulsar RXS J1708-40. 

trends exist between the pulse fractions and the brightness 
of a source, as we discuss below. 

For the case of two antipodal polar caps discussed here, 
the amplitude of pulsations from a neutron star of a given 
radius decreases with increasing polar-cap surface area and 
increasing mass. This is shown in Figure 4, where the max- 
imum pulse fraction, defined as 



PF 



F($)|max-F($)|r 
F^lmax + JX*)!, 



(16) 



which occurs for an orthogonal rotator (a = (3 = 90°), is 
plotted against the fractional surface area of the emitting 
region, for different neutron-star radii. 

As it has already been discussed extensively in the liter- 
ature, for the isotropic beaming we use here, general rel- 
ativistic light bending suppresses the pulse fraction below 
<~ 35%, even for infinitesimally small polar-cap sizes. As 
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Fig. 6. — Fraction N ^, s /Ntot of observed systems with pulse fraction larger than a given threshold as a function of the fractional size of 
the polar caps, for different neutron-star radii. 



a result, detection of a pulsation from a neutron star with 
a larger amplitude severely constrains any such models 
of emission from the stellar surface. However, these con- 
straints may become even tighter if the observed brightness 
of the neutron star is also taken into account. For the same 
local radiation spectrum, small polar-cap sizes correspond 
to large pulse fractions but weak radiation fluxes and vice 
versa. As a result, in this case, the maximum pulse frac- 
tion for a bright source is significantly smaller compared 
to the maximum pulse fraction for a faint source. 

Inferring the fractional surface area of the emitting re- 
gion, and thus the source brightness, requires an a priori 
knowledge of the orientation of the polar caps and the ob- 
server with respect to the rotation axis, which are almost 
always unknown. For example, a given observed source 
brightness may be the result of a small polar cap viewed 
from a favorable orientation (a ~ (3 ~ 0°) or of a larger 
polar cap viewed from an unfavorable orientation (a ~ 0° , 
(3 ~ 90°). Therefore, since the observed brightness does 
not have an one-to-one correspondence with the fractional 
emitting surface area, the constraints plotted in Figure 4 
cannot be compared directly with observations. However, 
although both configurations in the above example may 
produce the same pulse-average flux as measured by the 
observer, the first configuration typically produces a signif- 
icantly smaller pulse fraction. As a result, the combination 
of these two properties, namely the source brightness and 
pulse fraction, provide us with a useful diagnostic tool, as 
we discuss below. 

We can simultaneously account for all the above effects 
if we assume a particular model of the local radiation 
spectrum and search for the maximum pulse fraction as 
a function of the brightness of the source measured by an 
observer at infinity. As an example, we assume that the 
local radiation spectrum at each point on the polar caps 
is that of a blackbody of temperature Tns with isotropic 
beaming. An observer at infinity, would observe a radia- 
tion flux Foo, and making the assumption that the emis- 
sion is spherically symmetric, would infer a luminosity 

Lo. = A-kcPF^ . (17) 

According to equation 0] , this luminosity is 

Loo = ^R 2 NS aT* s (V=g^) 2 A(a, 0, 7 , r) , (18) 



where Tns = T^/ Soo)- We define an apparent bright- 
ness as 

/ tofFn \( ^ y 4 
~ ^10 36 ergs- 1 J ^0.5 keV/ ' 1 ' 

which is a function of only observed quantities and is in- 
dependent of the neutron-star temperature in our model 
calculations. We use this quantity as a measure of the 
brightness of the star at infinity and, in Figure 5, we plot 
against it the maximum pulse fraction for a 10 km neutron 
star emitting a blackbody spectrum. This figure shows 
clearly that, even though emission from a neutron-star 
surface can in principle produce pulse fractions as large 
as ~ 30% (for very small hot spots), it can never produce, 
e.g., a pulse fraction of ;> 10% simultaneously with an in- 
ferred luminosity of ;> 10 36 erg s -1 for a 0.5 keV blackbody 
temperature. 

4. DISCUSSION 

In this paper we have shown that the inferred proper- 
ties of thermally emitting neutron stars, i.e., their emitting 
surface areas and pulse fractions, are significantly affected 
by the three-dimensional geometry of the systems, phase- 
averaging effects, and general relativistic light bending. 
The uncertainties introduced by these effects can be sig- 
nificant (with a spread of <; 2), especially at the limit of 
small emitting surface areas. However, even when most of 
the neutron star surface is emitting, the uncertainties in 
estimating the neutron-star radius are of the same order 
as the ~ 5% level required (Lattimer & Prakash 2000) for 
constraining the equation of state of neutron-star matter. 

We have also argued that for a given local radiation 
spectrum emerging from a bright spot on the stellar sur- 
face, faint sources can give rise to larger pulse fractions 
than brighter sources. The maximum pulse fraction as a 
function of the brightness of the neutron star, measured 
by I (cf. eq. provides a diagnostic tool for distin- 

guishing between different emission models. These derived 
constraints depend both on the local radiation spectrum 
(e.g., thermal versus non-thermal emission) and its beam- 
ing and are, therefore, different for different classes of mod- 
els. Moreover, both quantities are measurable from spec- 
tral and timing observations and can be directly compared 
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to the calculated constraints. The high observed pulse 
fractions of AXPs (see, e.g., Chakrabarty et al. 2000) and 
X-ray bursters (Strohmayer et al. 1998) have been shown 
to strongly constrain emission models and such constraints 
can become only tighter when the apparent luminosities 
of these systems are also taken into account. As an ex- 
ample, we consider the source RXS J1708— 40, which has 
been identified as an AXP (Sugizaki et al. 1997). Fitting a 
blackbody spectrum to ASCA observations of this source 
gives a temperature of = 0.41 ± 0.03 keV and a black- 
body flux of Air^Foe = 3.16 x 10 35 erg s _1 (Sugizaki et 
al. 1997; Chakrabarty et al. 2000). The inferred apparent 
brightness of this source is, therefore, I = 0.8 and com- 
bined with the observed ~ 50% pulse fraction (Sugizaki 
et al. 1997) is inconsistent with isotropic thermal emission 
from the neutron-star surface (see Fig. 5). 

Finally, our results also have a number of important im- 
plications for soft X-ray surveys, in supernova remnants, 
for young, cooling neutron stars that emit thermally. For 
a given spectral shape, the brightest sources, which would 
be more easily detectable, correspond to smaller pulse frac- 
tions, while strong pulsations are only expected from dim- 
mer sources. This effect is shown in Figure 6, where the 
fraction of systems N b s /N t ot with a pulse fraction at in- 
finity larger than a threshold PFq is plotted against the 
fractional emitting surface area. For this purpose, we as- 
sume a random orientation of the magnetic inclination and 
the inclination to the observer for a sample of systems and 
define 

n v, f v/2 r^ 2 

-r^ 1 (PF )= / X[PF(a,/3);PF }smada sin/3d/3 , 

"tot Ja=0 J/3=Q 

(20) 



where the step function X(PF; PFo) is defined such that 
it is unity when PF > PF and zero at all other values of 
the pulse fraction. For realistic neutron star masses and 
radii (p ~ 2 — 3), only a very small fraction of sources 
shows pulsations that are detectable at a significant level 
(;> 10 — 20%) and this fraction drops rapidly with increas- 
ing apparent luminosity. 

If the central source in the remnant Cas A is a young, 
cooling neutron star, its surface brightness distribution 
cannot be uniform, as inferred from fitting thermal mod- 
els to the observed countrate spectra (Pavlov et al. 2000; 
Chakrabarty et al. 2000). However, as Figure 6 shows, 
this property is not inconsistent with the ~ 30% upper 
limit on its pulse fraction (Chakrabarty et al. 2000). For 
the assumed isotropic beaming of the emerging radiation 
and the very small fractional surface area inferred for the 
central source in Cas A (Pavlov et al. 2000; Chakrabarty 
et al. 2000), less than half of the systems would have been 
detected with pulse fractions higher than the detection 
threshold, if p = 4. For more realistic neutron-star prop- 
erties (p = 2, 3), no system would show such a large pulse 
fraction. Therefore, the absence of detectable pulsations 
from this source is not a strong argument against its iden- 
tification with a spinning neutron star. 

We thank Deepto Chakrabarty, Lars Hernquist, and 
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knowledges the support of a post-doctoral fellowship from 
the Smithsonian Insitution. F. O. acknowledges support of 
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